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Abstract 

We examine the interference 7 X box corrections to parity-violating elastic electron-proton scattering in the 
light of the recent observation of quark-hadron duality in parity-violating deep-inelastic scattering from the 
deuteron, and the approximate isospin independence of duality in the electromagnetic nucleon structure 
functions down to « 1 GeV^. Assuming that a similar behavior also holds for the jZ proton structure 
functions, we find that duality constrains the jZ box correction to the proton’s weak charge to be 3?e = 

(5.4 ± 0.4) X 10“^ at the kinematics of the Qweak experiment. Within the same model we also provide 
estimates of the fZ corrections for future parity-violating experiments, such as MOLLER at Jefferson Lab 
and MESA at Mainz. 
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1. Introduction 

Parity-violating precision measurements have for many years provided crucial low-energy tests of the 
Standard Model. Early efforts such as the E122 experiment at SLAG m firmly established the SU(2)xU(l) 
model as the theory of the unified electroweak interactions. Modern-day experiments use parity violation 
to probe physics beyond the Standard Model. One of the most recent parity-violating measurements is the 
Qweak experiment at Jefferson Lab [3], which aims to measure the proton’s weak charge to 4% accuracy. With 
an initial analysis of a subset of the data already reported [ 1 ] , the analysis of the full data set is expected in 
the near future. 

For the precision requirements of the Qweak experiment, the weak charge of the proton, defined at tree 
level as = 1 —4sin^ 9w, must also include radiative corrections. Including these corrections at the 1-loop 
level, the weak charge can be written as [5] 

Qw — (1 + Ap-I-Ae) — 4sin^ 0 u/(O)-I--I-Dwiv + Glzz + □ 7 ^( 0 ), ( 1 ) 

where sin^ 6^ ( 0 ) is the weak mixing angle at zero momentum transfer, and the electroweak vertex and 
neutral current correction terms Ap, Ag and Ag have been calculated to the necessary levels of precision [5]. 
The weak box corrections □ww are dominated by short-distance effects and can also be computed 

perturbatively to the required accuracy. 

On the other hand, the final term in Eq. Q, the 7 X box contribution, depends on both short- and 
long-distance physics and therefore requires nonperturbative input. Gonsiderable attention has been given to 
the analysis of this term, for both the vector electron-axial vector hadron coupling to the Z, (which is 
relevant for atomic parity violation experiments) [SI [7l [HI [9] , and the axial electron-vector hadron coupling, 
'^Xz (which because of its strong energy dependence makes important contributions to the Qweak experiment) 
nni dn ini (mill. The most accurate technique to evaluate the latter is a dispersion relation. While 
constraints from parton distribution functions (PDFs) and recent parity-violating deep-inelastic scattering 
(PVDIS) data [T5||T6| provide a systematic way of reducing the errors on this correction [T4j, some uncertainty 
remains about the model dependence of the low-Q^ input. 
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The E08-011 electron-deuteron PVDIS experiment at Jefferson Lab not only allowed an accurate deter¬ 
mination of the C 2 q electron-quark effective weak couplings |16j . but also presented the first direct evidence 
for quark-hadron duality in 'yZ interference structure functions, which was verified at the (10-15)% level for 

down to « 1 GeV^ [15j . In general, quark-hadron duality refers to the similarity of low-energy hadronic 
cross sections, averaged over resonances, with asymptotic cross sections, calculated at the parton level and 
extrapolated to the resonance region. It is manifested in many different hadronic observables m and was 
first observed in deep-inelastic scattering (DIS) by Bloom and Gilman [TSl [TO]. Subsequent studies have 
quantified the validity of duality for various spin-averaged and spin-dependent electromagnetic structure 
functions, as well as in neutrino scattering and for different targets establishing 

the phenomenon as a general feature of the strong interaction. 

Furthermore, recent analysis of moments of the free neutron electromagnetic structure function |28j has 
demonstrated that duality in the lowest three neutron moments is violated at a similar level 10%) as in the 
proton for > 1 GeV^ (TO] [^ |2S| . This suggests that the isospin dependence of duality and its violation 
is relatively weak. It is reasonable therefore to expect that duality may also hold to a similar degree for the 
jZ structure functions, which are related to the electromagnetic structure functions by isospin rotations. 

In this paper we discuss the extent to which quark-hadron duality in yZ structure functions can provide 
additional constraints on the corrections, and in particular the contributions from low hadronic final 
state masses W and ~ 1 GeV^. In Sec. we illustrate the realization of duality in the moments of 
the proton and neutron electromagnetic structure functions using empirical parametrizations of data in the 
resonance and DIS regions down to = 1 GeV^. Motivated by the approximate isospin independence of 
duality in electromagnetic scattering from the nucleon, in Sec. 1^ we explore the consequences of duality in 
the yZ structure functions for the energy dependence of the correction, and especially the limits on 
its overall uncertainty. Finally, in Sec. we summarize our findings and discuss their implications for the 
analysis of the Qweak experiment as well as future parity-violating experiments such as MOLLFR at Jefferson 
Lab [TO] and MFSA at Mainz [30] . 


2. Duality in electromagnetic structure functions 


Historically, the observation of duality in inclusive electron scattering [TOl [TO] predates the development 
of QCD and was initially formulated in the language of finite-energy sum rules. Within QCD, duality was 
reinterpreted within the operator product expansion through moments of structure functions m, with duality 
violations associated with matrix elements of higher twist (HT) operators describing multi-parton physics. 
The extent to which inclusive lepton-nucleon cross sections can be described by incoherent scattering from 
individual partons through leading twist (LT) PDFs can be quantified by studying the dependence of 
the structure function moments. At low corrections to the LT results arise not only from multi-parton 
processes, but also from kinematical target mass corrections (TMGs), which, although 1/Q^ suppressed, arise 
from LT operators. To isolate the genuine duality-violating HT effects, one can consider Nachtmann moments 
of structure functions [32] , which are constructed to explicitly remove the effects of higher spin operators and 
the resulting TMGs. 

Specifically, the Nachtmann moments of the Fi and F 2 structure functions are defined as [331 [TO] 
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is the Nachtmann scaling variable [TO] [TO], with x = -|- Q^) the Bjorken scaling variable, 

= 1 -I- AM^x^/Q^, and M the nucleon mass. The variable pn is given by 
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Figure 1: The proton (left panels) and neutron (right panels) electromagnetic (top) s,nd (bottom) structure function 
moments. The total Nachtmann moments (black solid lines) include contributions from the resonance {W'^ ^ 6 GeV^, blue 
dot-dashed lines) and DIS > 6 GeV^, green dotted lines) regions, as well as the elastic contributions (gray dashed lines), 

and are compared with the Cornwall-Norton moments of the LT structure functions (red long-dashed lines). 


and vanishes in the —)■ oo limit. In that limit the moments approach the standard Cornwall-Norton 

moments [55] . 

= f dxx^-^F,{x,Q^), i = l,2. (6) 

Jo 

At finite while the moments depend only on the F 2 structure function, the moments have 
contributions from both the Fi and F 2 structure functions. Because the latter contribution is proportional 
to rin, it vanishes at large so that the moments are generally dominated by the Fi structure function 
at large Q^. 

Duality in unpolarized electron-nucleon scattering has been studied most extensively for the electromag¬ 
netic F 2 structure function [201 [HIUS], and to a lesser extent for the Fi (or longitudinal Fl) structure function 
iniET]. The latter is generally more difficult to access experimentally, as it requires precise longitudinal- 
transverse separated cross section measurements, or equivalently the (Jl!'Xt cross section ratio. In Fig. [^the 
workings of duality in the n = 2 Nachtmann moments of the proton and neutron Fj'^ and F'^'^ structure 
functions are illustrated over the range 1 ^ ^ 8 GeV^. For the low-kF^ contributions, IF^ ^ 6 GeV^, 

the resonance-based fit to the electromagnetic structure function data from Christy and Bosted [55] is used. 
For the DIS region at higher values, > 6 GeV^, this is supplemented by the ABM global QCD fit 
[39] to high-energy data, which includes LT, TMC and HT contributions. Since LT evolution is logarithmic 
in Q^, at large the moments are predicted to become flat in InQ^. While the individual resonance and 
DIS region contributions, as well as the elastic (W = M) component, are strongly dependent in the 
region of low shown in Fig. remarkably their sum exhibits only very mild dependence down to 
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« 1 GeV^. This is the classic manifestation of duality observed by Bloom and Gilman nsms], in which 
the total empirical moments resemble the LT contributions down to surprisingly low momentum scales. Note 
that since the Nachtmann moments are constructed to remove higher spin operators that are responsible for 
TMGs, in the absence of HTs one would expect the Nachtmann moments of the total structure functions to 
equal the Gornwall-Norton moments of the LT functions, /i^"'^(LT + TMC) = Mj*'"^(LT) [ID] . 

This expectation is clearly borne out in Fig. where the total and moments are very similar 
to the moments computed from the LT PDFs. For the proton structure functions, the average violation of 
duality in the range 1 ^ ^ 2.5 GeV^ is 3% and 4% for the and structure functions, respectively, 

with the maximum violation being « 5% and « 10% at the lower end of the range. For the neutron the 
maximum violation is slightly larger, with the LT F^"^ and F^^ moments being « 14% and « 10% smaller 
than the full results, although the average over this range is 5% and 8%, respectively. This is consistent 
with several previous phenomenological analyses of high-energy scattering data which have found 

no indication of strong isospin dependence of HT corrections. Following Ref. [5S|, we assign a 5% error on 
the proton and F^'^ structure functions, and a larger, 10% error on the neutron structure function [il] . 
reflecting the additional nuclear model dependence in extracting the latter from deuterium data |45j . For the 
elastic contribution a 5% uncertainty is assumed for the total elastic structure functions from Ref. [46]. For 
higher moments (n > 2), which are progressively more sensitive to the high-a; (or low-VF) region, the degree 
to which duality is satisfied diminishes at lower values m- 


3. Duality in jZ structure functions and implications for 


In contrast to the electromagnetic structure functions which have been studied extensively for many years, 
experimental information on the interference "fZ structure functions is for the most part nonexistent. Some 
measurements of F^^ and xF^^ have been made at very high at HERA [IS], where the yZ contribution 
becomes comparable to the purely electromagnetic component of the neutral current. However, no direct 
measurements of F^^ and F^^ for the proton exist in the ^ few GeV^ range relevant for the evaluation 
of the "fZ box correction to dH . 

In principle, the computation of the imaginary part of the correction to the proton’s weak charge at 
a given incident energy E requires knowledge of the ^Z structure functions over all kinematics. 
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where a is the running electromagnetic coupling evaluated at the scale and Mz is the Z boson mass. 
The W'^ range covered in the integral lies between the inelastic threshold, = (M -|- and the total 
electron-proton center of mass energy squared, s = -|- 2ME, while the integration range is from 0 

up to Qmax = 2MF(1 — VF^/s). (The small mass of the electron is neglected throughout.) The real part of 
the jZ box correction which enters in Eq. (IR can then be determined from the imaginary part through an 
unsubtracted dispersion relation [inunKniiiiiii, 

n p nOO 1 

5ien^2(F) = —Vj^ (8) 

where V is the Gauchy principal value integral. While the dispersion relation Q is valid only for forward 
scattering, because the Qweak experiment is performed at a small scattering angle « 6°, in practice it provides 
a very good approximation. 

Note that at high and large E, the total correction can also be expressed in terms of the 

moments of the F^^ and F^^ structure functions by switching the order of the integrations in Eqs. ^ and 
([^ and expanding the integrand in powers oi x^/Q^ |5]. The higher order terms in 1/Q^ are then given 
in terms of higher moments of the structure functions. The expansion in Ref. jS] was performed in terms 
of the Cornwall-Norton moments, but the expansion could also be generalized to the Nachtmann moments 
in Eqs. and @. However, because this approximation neglects contributions from the low-lF region, it 
is appropriate only for DIS kinematics and is not directly applicable for the present application, where the 
integrals are dominated by contributions at low and IF^. In particular, as we discuss below, at energy 
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Figure 2: Kinematic regions contributing to the integrals in the AJM model. Region I (blue) includes the nucleon resonance 
region at low and Q^\ Region II (red) encompasses the low-Q^, high-IK^ region described by Regge theory; and Region III 
(green) is the deep-inelastic region characterized by LT PDFs. The shaded band between = 1 and 2.5 GeV^ represents the 
extension of Region III from its previous boundary in Ref. |14| (Q^ = 2.5 GeV^) to its current reach {Q^ = 1 GeV^). 


E ^ 1 GeV, approximately 2/3 of the integral comes from the traditional resonance region W < 2 GeV and 

< 1 GeV^. In contrast, the contribution from the DIS region for ly > 2 GeV and > 1 GeV^ is « 13% 
at this energy. 

In Refs. [miSn] the and structure functions were computed from the phenomenological Adelaide- 
Jefferson Lab-Manitoba (AJM) parametrization. This is based on the electromagnetic structure functions 
described in Sec. but appropriately rotated to the 'yZ case according to the specific W'^ and region 
considered, with the rotation parameters constrained by phenomenological PDFs [14] and recent PVDIS data 
[T5l ITH] . In the AJM model the integrals over and in Eq. ([^ are split into three distinct regions, 
characterized by different physical mechanisms underlying the scattering process. In each region the most 
accurate parametrizations or models of and F^^ available for the appropriate kinematics are used. 

In the present analysis, we define the W'^ and regions as illustrated in Fig. “Region I” (low Q^, low 
encompasses 0 ^ < 10 GeV^ for < 4 GeV^, and 0 ^ < 1 GeV^ for 4 < < 9 GeV^, 

using the 77 —>■ jZ rotated Christy-Bosted parametrization [38] of the resonance + background structure 
functions. For “Region 11” (low Q^, high W'^), the vector meson dominance + Regge model of Alwall and 
Ingelman [JS] is used over the range 0 ^ ^ 1 GeV^ and > 9 GeV^. Finally, for “Region III” (high 

high W^) the perturbative QGD-based global fit from Alekhin et al. (ABM) [ 35 ] is used for > \ GeV^ 
and > 4 GeV^, which includes LT as well as subleading 1/Q^ TMG and HT contributions. For x = 1, 
the elastic contributions to the structure functions are computed using the form factor parametrizations from 
Ref. [46] . 

While the uncertainties on the 7 Z structure functions in Region III are small — typically a few %, 
reflecting the errors on the PDFs from which they are constructed through the simple replacement of quark 
charges —>■ — the uncertainties in F^^ and F^^ are expected to be larger at lower and Q^. In the 
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previous analyses of the 7 Z correction misn], the PDF-based description was limited to > 2.5 GeV^ (and 
> 4 GeV^). Motivated by the observation of duality in the proton and neutron and structure 
functions, and in PVDIS from the deuteron, as discussed in Sec. we further assume the approximate 
validity of duality in the jZ proton structure functions and extend the QCD description of Region III down 
to = 1 GeV^. Lowering the boundary of the DIS region, which is well constrained by leading twist PDFs, 
to smaller decreases the contribution from Regions I and II, and hence reduces the model uncertainty on 
the 77 —>■ rotation of the structure functions in this region. 

Within the AJM 7 Z structure function parametrization, the most uncertain elements are the con¬ 
tinuum parameters used to relate the high-mass, non-resonant continuum part of the ^Z transverse and 
longitudinal cross sections to the 77 cross sections in the generalized vector meson dominance model [13 (51]. 
The Kq parameters are htted by matching the 7 Z to 77 cross section ratios with the LT structure function 
ratios at = 1 GeV^, 


\i^c) 
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where the longitudinal structure function is related to the F\ and F 2 structure functions by = 
P^F 2 — 2xFi [H]. (Note that, consistent with the duality hypothesis, we use the LT structure functions 
in Region III rather than the total structure functions that may include the small subleading contributions 
[33 •) The resulting fit values, 

= 0.36 ±0.15, = 1.5 ±3.1, (10) 

are obtained by averaging over the parameter determined from 10 fits with the ratios in Eq. ^ matched 
at between = 4 GeV^ and 13 GeV^. These values are then used to compute the ^Z structure functions 
in the dispersion integral for 1 ^ ^ 10 GeV^ and ^ IF^ ^ 4 GeV^. To allow for stronger violations 

of duality at lower Q^, the uncertainties on are inflated to 100 % for the region 0 ^ < 1 GeV^ for all 

IF^. In the numerical calculations the uncertainties on the proton jZ structure function parametrizations 
are taken to be the same as those used in the calculation in Ref. [14] , and a 5% uncertainty is assumed 
for the nucleon elastic contributions. 

Using the jZ structure functions obtained from the newly fitted values, the correction is 

displayed in Fig. [^ as a function of beam energy, with a breakdown of the individual contributions from 
different regions given in Table[^ At the incident beam energy F = 1.165 GeV of the Qweak experiment, the 
total correction is found to be 


?fieDT = (5.4±0.4) X 10"^ (11) 

This is in good agreement with the value = (5.57 ± 0.36) x 10“^ found in the previous analysis [14]. 

In particular, even though the values of the continuum rotation parameters in the earlier fit were somewhat 
different {kq = 0.65 ± 0.14 and Kq = —1.3 ± 1.7 with matching to the total DIS structure functions at 
= 2.5 GeV^), the central value of remains relatively unaffected. 


Table 1: Contributions to from Regions I, II and III, and the total, at the kinematics of the Qweak? MOLLER, and 

MESA experiments. 




^eUXz (xlO-3) 



Qweak 

MOLLER 

MESA 

Region 

(F =1.165 GeV) 

(F = 11 GeV) 

(F = 0.18 GeV) 

I 

4.3 ±0.4 

2.5 ±0.3 

1 . 0 ± 0.1 

II 

0.4 ±0.05 

3.2 ±0.5 

0.06 ± 0.01 

III 

0.7 ±0.04 

5.5 ±0.3 

0.1 ± 0.01 

Total 

5.4 ±0.4 

11.2 ±0.7 

1.2 ± 0.1 
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Figure 3: Energy dependence of the 7 Z box correction, Qw‘ contributions from various regions in and 

(Regions I, II and III) are shown separately, as is the total (solid curve). The dashed vertical lines indicate the beam energies 
of the various parity-violating experiments {E = 0.18 GeV for MESA | 30| . E = 1.165 GeV for Qweak [I]) E = 11 GeV for 
MOLLER [29]. 


The largest contribution to 3fie at the Qweak energy is still from Region I, which makes up « 80% of 
the total, with its error dominating the total uncertainty. Of this, « 2/3 is from the traditional resonance 
region W'^ < 4 GeV^ (of which 61% is from < 1 GeV^ and 6 % from > 1 GeV^), and « 13% is 
from < 1 GeV^ and 4 < < 9 GeV^. The contributions from Regions II and III are « 7% and 

« 13%, respectively, of the total at the Qweak energy, but become more important with increasing energy. 
Interestingly, the modified boundary for Region III results in a somewhat smaller contribution from 
Region II (0.4 x 10“^ compared with 0.6 x 10“^), while the Region III contribution has doubled (0.7 x 10“^ 
compared with 0.35 x 10“^) relative to that in Ref. m- In effect, moving the boundary from 2.5 GeV^ 
to 1 GeV^ shifts « 6 % of the total correction 3?e from Regions I and II to Region III. 

Furthermore, since the structure functions at < 1 GeV^ depend on because the k values are 
refitted at = 1 GeV^, duality also indirectly affects the low-Q^ contribution. Therefore, although duality 
is formally used only down to = I GeV^, the constraint influences the jZ calculation below 1 GeV^ as 
well, as the matching now is to a more reliable ■jZ cross section at that point. 

While we have assumed the validity of duality for the and structure functions down to = 
I GeV^, the possible violations of duality have a minor effect on the analysis. Even if one takes the maximum 
violation of duality (« 14%) in the 77 structure functions seen in Fig. at the lowest over the entire 
1 ^ ^ 2.5 GeV^ range, the error introduced into the total from duality violation is < 0.1%. 

Overall, compared with Ref. [T3] the total relative uncertainty increases marginally, from 6.5% to 7.4%, 
despite the rather more conservative estimates of the structure function uncertainty for ^ 1 GeV^ through 
the inflated errors on . Note that the same 100% uncertainties are used in the transformation of the 
vector meson dominance model gsun] in Region II. For Region III, the LT F^^ and F^^ structure functions 
are assigned a 5% uncertainty for ^ 2.5 GeV^, which is increased linearly to 10% at = 1.0 GeV^. 

Since the electromagnetic structure functions are reasonably well approximated by the LT results even 
below the traditional resonance-DIS boundary of = 4 GeV^, we also examine the effect of lowering the 
cut into the peripheral resonance region down to = 3 GeV^. In this case the contribution from 
Region III increases to 0.9 x 10“^, while that from Region I correspondingly decreases to 4.2 x 10“^, hence 
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leaving the total essentially unchanged. 

At the higher E = 11 GeV energy of the planned MOLLER experiment at Jefferson Lab [23], the DIS 
region contributes about half of the total, = (11.2 ±0.7) x 10“^, with Regions 1 and II making up the 

other 50%. This again agrees well with the earlier determination = (11.5±0.8) x 10“^ from Ref. [50] . 

On the other hand, for the possible future MESA experiment in Mainz m at a lower energy, if = 0.18 GeV, 
the bulk of the contribution still comes from Region I, but is reduced by a factor of ~ 4 compared with the 
correction at the Qweak energy. 

4. Conclusion 

Quark-hadron duality is one of the most remarkable phenomena ever observed in hadronic physics. While 
some aspects of global duality can be formulated in the language of QGD, such as the relation between the 
scale independence of structure function moments and the size of higher twists, the detailed workings of local 
duality, for specific regions of W'^ or x, are not well understood from first principles. Nevertheless, there are 
many marvellous practical applications to which duality can be put. For example, the high-energy behavior 
of hadronic cross sections can be used to predict averages of resonance properties; and, conversely, low-kE^ 
data, suitably averaged, can be utilized to constrain LT parton distributions in difficult to access kinematic 
regions. 

The latter category appears the most promising approach at present, with several global PDF analyses 
|39||45l|43| extending their coverage down to lower > 1 GeV^) and {W"^ > 3 GeV^) values than 

in traditional LT analyses. This not only increases considerably the available data base for PDF fitting, it is 
also one of the few ways currently available to study PDFs at high a; ^ 1. 

The main implication of duality for the current analysis is the extension of the LT description of jZ 
structure functions to lower = 1 GeV^, than in previous work [14]. This serves to reduce the size 

of the contribution from Region I, which has the largest uncertainty associated with the behavior of the 
yZ strncture functions at low and W^. To account for the possible model dependence of the 77 —>■ jZ 
structure function rotation and the violation of duality at low , we have assigned rather conservative errors 
on and in this region. This is reflected in the increased uncertainty on this contribution compared 
with our previous analysis which is somewhat offset by the larger contribution from Region III that is 
well constrained by PDFs. 

The final result of = (5.4 ± 0.4) x 10“^ is consistent with Ref. m, but with a slightly larger 

relative uncertainty, which comes almost entirely from Region 1. It also agrees with the central value from 
Ref. [13j . although the error there is « 5 times larger, which in view of our current analysis appears to 
be somewhat overestimated. Our findings suggest that with the constraints from existing PVDIS data and 
PDFs, and now with the further support from quark-hadron duality, the overall uncertainty in the estimate 
of the "fZ box correction is well within the range needed for an unambiguous extraction of the weak charge 
from the Qweak experiment. 

Further reduction of the uncertainty on the "fZ correction will come from new measurements of PVDIS 
asymmetries on the proton, particularly at the low and values that are most relevant at the 
Qweak energy. These will also be useful in constraining the ■jZ contribution at the much lower energy 
if = 0.18 GeV of the MESA experiment [30], where we find the correction to be » 4 times smaller but even 
more dominated by Region 1. In contrast, for the MOLLER experiment at the higher E = 11 GeV energy 
the dispersion integral is dominated by the DIS region, which although contributing to a larger overall 
correction, is better determined in terms of PDFs. These new experiments hold the promise of allowing the 
most precise low-energy determination of the weak mixing angle to date, and providing a unique window on 
possible new physics beyond the Standard Model. 
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